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1.  INTRODUCTION 

This  study  is  concerned  with  understanding  the  nonlinear  coupling 
which  occurs  between  the  fluid  dynamics,  kinetics  and  the  optical  reso¬ 
nator  of  a  cw  chemical  laser.  The  motivation  for  this  work  is  the  neces¬ 
sity  of  predicting  the  power  and  power  spectral  performance  of  the  laser. 
The  significant  effects  of  rotational  nonequilibrium  on  the  performance  of 

chemical  lasers  have  been  described  in  the  HF  cw  case  by  Sentman^  and  in 

2 

the  pulsed  case  by  Hough  and  Kerber  .  The  DF  cw  case  has  been  treated  by 
3  _  _ 

Hall  and  the  HF  single  line  amplifier  by  Skifstad  .  The  agreement  be¬ 
tween  these  calculations  and  the  experimental  data  showed  that,  when  ro¬ 
tational  nonequilibrium  is  included  in  the  model,  several  J  transitions 
within  a  given  vibrational  band  lase  simultaneously,  most  of  the  power  is 
contained  in  those  J  transitions  which  are  preferentially  populated  by  the 
pumping  reaction,  and  the  power  is  decreased  by  20-30%  from  the  rotational 
equilibrium  value.  A.11  of  these  studies  employed  a  Fabry-Perot  model  for 
the  optical  cavity,  whereas  practical  lasers  usually  extract  power  with 
an  unstable  resonator.  The  effect  of  the  optical  model  on  cw  laser  per¬ 
formance  was  examined  by  Sentman"’ .  Calculations  were  performed  for  both 
rotational  equilibrium  and  nonequilibrium  cases  for  a  confocal  unstable 
resonator  with  50%  geometric  outcoupling.  For  the  rotational  equilibrium 
case,  the  solution  resulted  in  most  of  the  power  being  contained  in  one 
line.  For  the  rotational  nonequilibrium  case,  the  power  is  distributed 
over  eleven  lines,  and  the  Fabry-Perot  and  unstable  resonator  power  spec¬ 
tral  distributions  are  almost  identical,  the  major  difference  being  that 
all  lines  whose  saturated  gain  region  does  not  extend  into  the  central 
portion  of  the  unstable  resonator  are  cut  off.  For  the  rotational  non¬ 
equilibrium  case,  fluctuations  in  cw  power  were  shown  to  occur  on  lines 
whose  saturated  gain  zone  does  not  fill  the  unstable  resonator.  The 


amplitude  of  the  fluctuation  was  determined  by  the  fraction  of  the  resonator 
filled  by  the  saturated  gain  zone  of  the  oscillating  line,  and  the  frequency 
of  the  fluctuation  was  determined  by  the  location  of  the  intensity  peak 
of  the  oscillating  line. 

These  studies  have  demonstrated  that: 

1.  Rotational  nonequilibrium  plays  a  major  role  in  determining  the 
power  spectral  performance  of  the  laser.  This  phenomenon  must 
be  included  in  any  model,  Fabry-Perot  or  unstable  resonator, 
which  is  to  predict  the  power  spectral  performance  of  the  laser. 

2.  Rotational  equilibrium  models  over  predict  the  total  power  and 
maximum  intensity. 

3.  For  rotational  nonequilibrium,  the  Fabry-Perot  and  unstable  reso¬ 
nator  power  spectral  distributions  are  identical. 

4.  Rotational  nonequilibrium  is  the  mechanism  responsible  for  a  non¬ 
linear  coupling  which  may  occur  between  the  saturated  gain  and  the 
resonator  geometry  resulting  in  large  amplitude  fluctuations  on 

a  given  line. 

Since  the  computation  time  for  a  computer  model  of  a  chemical  laser 
2 

depends  upon  N  ,  where  N  is  the  number  of  dependent  variables,  the  neces¬ 
sity  of  including  rotational  nonequilibrium  in  the  model  will  result  in 
a  significant  increase  in  run  time  because  the  number  of  dependent  vari¬ 
ables  will  increase  from  about  30  to  80  or  more.  Since  most  of  the  pre¬ 
ceding  studies  were  carried  out  using  very  simplified  kinetic,  fluid  dy¬ 
namic  models,  these  results  are  qualitative  rather  than  quantitative. 

Thus,  the  main  thrust  of  the  present  work  was  to  modify  the  qualitative, 
but  computationally  efficient,  rotational  nonequilibrium  cw  chemical  laser 
model  developed  by  Sentman^’^  so  that  it  would  give  quantitative  predic¬ 
tions  of  the  power  spectral  performance  of  the  laser. 


_ 
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Since  coupled  kinetic-physical  optics  calculations  are  iterative  in 
nature  and  computationally  expensive  300  sec/ i teration  on  the  CYBLR  175 
with  Sentman's  model^  with  convergence  occurring  in  about  15  iterations), 
a  two  pronged  approach  was  employed.  First,  the  modification  of  the  kinetic- 
fluid  dynamic  model  was  carried  out  using  a  Fabry-Perot  model  *  for  the  reso¬ 
nator  (a  typical  case  with  Sentman's  model1  takes  30  sec  on  the  CYBFR  175). 
Simultaneously,  a  physical  optics  model  supplied  by  AFWL  was  coupled  to 
Sentman's  rotational  nonequilibrium  model1  and  several  of  the  previous  cases 
run  to  verify  those  results.  The  AFWL  optics  model  was  then  coupled  to 
the  new  kinetics-fluid  dynamic  model  which  has  been  developed.  In  Section 
II,  the  kinetic-fluid  dynamic  model  which  has  evolved  is  described;  the 
comparison  with  the  results  from  a  sophisticated  chemical  laser  model,  the 
Bell  Aerospace  Textron  Blaze  II  model0,  for  two  different  cases  is  given 
in  Section  III.  In  Section  IV,  the  comparison  of  the  AFWL  physical  optics 
model  with  the  original  physical  optics  model S  are  presented  and  the  status 
of  the  coupled  AFWL  physical  optics  and  the  new  kinetics  is  described. 

Several  conclusions  and  recommendations  for  future  studies  are  discussed 


11.  KlNli'l  ICS-1 LIJ Il»  DYNAMICS  MODliL 


The  main  thrust  of  this  work  was  to  modify  the  qualitative,  but  com¬ 
putationally  efficient,  rotational  nonequi l ib.  ium  cw  chemical  laser  model 
developed  by  Sentman1,5  so  that  it  would  give  quantitative  predictions  of 
the  power  spectral  performance  of  the  laser,  lor  computational  efficiency, 
the  number  of  dependent  variables  must  be  kept  to  a  minimum.  Based  on  ex¬ 
perience  with  that  model,  the  approach  used  in  this  study  was  to  include 
the  essential  chemical  kinetics  -  the  pumping  reactions,  major  deactivation 
reactions  and  rotational  relaxation  reactions  -,  and  to  include  the  fluid 
dynamic  effects  of  mixing  by  inputing  the  pressure,  temperature,  velocity, 
mass  flow  rate  remaining  in  the  primary,  mass  flow  rate  remaining  in  the 
secondary  and  the  length  of  the  active  media  as  functions  of  x  which  are  ob¬ 
tained  from  one  of  the  available  quasi-2D,  2D  or  3D  fluid  dynamic,  rota- 

6  7  8 

tional  equilibrium,  chemical  kinetic  laser  models  ’  ’  .  This  approach  was 
motivated  by  the  fact  that  changes  in  chemical  kinetics,  which  significantly 
affect  the  power  spectral  performance,  in  most  cases,  have  only  a  minor 
effect  on  the  flow  field  in  which  the  kinetics  takes  place.  Thus,  the  fluid 
dynamics  provides  the  environment  in  which  the  kinetics  occurs. 

Since  this  model  will  reproduce  the  rotational  equilibrium  results  when 
it  is  run  using  a  value  of  the  rotational  relaxation  rate  constant  which  is 
105  times  its  measured  value,  the  model  has  evolved  from  detailed  compari¬ 
sons  with  the  Blaze  II  model^  for  two  different  cases:  an  Aerospace  arc 

g 

driven  HF  power  spectral  experiment  and  the  Bell  Aerospace  Textron  CL  XI 
(15,30)  lasing  lines  experiment.  These  two  cases  were  chosen  because,  in 
addition  to  providing  power  spectral  data, and  in  the  CL  XI  case,  intensity 
distributions  on  each  line,  they  represent  two  extremes  of  fluid  dynamic 
mixing;  the  Aerospace  case  is  a  slow  mixer  whereas  the  CL  XI  is  a  relatively 
fast  mixing  nozzle.  This  should  ensure  that  the  polynomial  fits  used  for 
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the  fluid  variables  should  be  capable  of  handling  any  profiles  likely  to 
be  encountered. 


2.1  THE  RADIATION  FIELD 

For  purposes  of  developing  the  kinetic-fluid  dynamic  model ,  the  opti¬ 
cal  cavity  was  modeled  as  a  Fabry-Perot  resonator.  The  physical  optics 
model  for  an  unstable  resonator  is  described  in  Chapter  IV. 

When  rotational  nonequilibrium  is  included  in  the  model,  several  J 
transitions  within  each  vibrational  band  lase  simultaneously  with  the 
dominant  ones  being  determined  by  the  J  dependence  of  the  pumping  re¬ 
action,  Ref.  1.  Because  of  this,  the  treatment  of  rotational  nonequili¬ 
brium  effects  is  facilitated  by  using  a  more  exact  treatment  of  the  radia¬ 
tion  field  than  that  afforded  by  the  constant  gain  approximation.  The  de¬ 
scription  of  the  radiation  field  employed  is  that  developed  by  Sentman , 
Ref.  1  and  summarized  below. 

The  geometry  is  shown  in  Fig.  1,  where  plane  parallel  mirrors  form  a 
Fabry-Perot  cavity.  The  gas  between  the  active  medium  and  the  mirrors  is 
assumed  to  be  nonabsorbing  and  hence  the  mirror  boundary  conditions  can  be 
applied  at  Z  =  0  and  L.  Tie  radiative  transfer  equations  are 


+  V 


(2.1-1) 


-  a  (I 

v  v 


V 


(2.1-2) 


where  is  the  gain  coefficient  and  is  the  spontaneous  emission  term 
given  by 

Sv  =  2h2v4NApBLUnunULgL/c2guav  (2.1-3) 


nM  and  n.  are  the  number  of  moles  in  the  upper  and  lower  states  per  gram 
u  ^ 
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of  fluid,  p  is  the  density  of  the  fluid,  nUL  is  a  line  shape  function  that 
is  normalized  according  to 


nUL  dV  =  1 


(2.1-4) 


is  the  Einstein  coefficient,  g.  is  the  statistical  weight  of  state  i. 


LU  '  “i 

N  is  Avogadro's  number  and  h  is  Planck's  constant.  The  general  solution 


of  Eqs.  (2.1-1)  and  (2.1-2)  averaged  over  the  narrow  spike  at  Vq  %  v  }  is 


I  (Z)  =  I  (0)  exp (a  Z)  +  K  [exp (a  Z)  -  1] 
UL  UL 


I"(Z)  =  I*(E)  exp [a  (L-Z)] +  K  [exp(a  [L-Z])  -  1 
UL  UL 


(2.1-5) 

(2.1-6) 


where  the  average  of  over  v  is  approximated  by.  Ref.  1, 


K  = 


S  dv  < 

S  dv 

J  v 

J  pn 

a'T 


and 


\l  =  ^UL^UL^LU^L7^  nU  '  nL] 


(2.1-7) 


(2.1-8) 


O'  is  the  Stefan-Boltzmann  constant,  g.  =  2J.  +  1  and  the  flow  and 
chemistry  are  assumed  independent  of  Z. 

The  quantities  I+(0)  and  I~(L),  which  are  determined  by  requiring 
Eqs.  (2.1-5)  and  (2.1-6)  to  satisfy  the  mirror  boundary  conditions 
I+(0)  =  rQI-(0)  and  I"  (L)  =  rLI+(L),  are 

I+(0)  =  rQK[rLexp(2aL)  +  (l-rL)exp(aL)  -  1]  / [  1-r^exp  (2aL)  ] 

(2.1-9) 

I"(L)  =  r  K[r  exp(2aL)  +  (1-r  ) exp(aL) - 1] / [ 1 -r  r  exp(2aL)] 

L  O  O  O  L 


where  the  subscript  has  been  suppressed.  Witfi  Eqs.  (2.1-9)  and 


(2.1-10) 
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(2.1-10),  the  intensity  averaged  over  Z  is  given  by 


I  =  I+  +  I"  =  K[  (rQ  +  -  2rQr^[  1-ixL]  )exp(2aL) 


+  2(1  +  r  -  ro  -  r^)exp(aL)  -  2(1  +  aL) 


rQ  +  rLJ/{aL[l  -  ror^exp(2aL) ] } 


(2.1-11) 


With  this  expression  for  I,  the  Xraj  term  in  the  species  equations  for 
n^j  and  n^  is  given  by 


X  ,  =  a  I  /hN.v... 

radUL  VUL  VUL  A  UL 


12.1-12) 


The  power  contained  in  the  transition  UL  is  obtained  by  integrating 


dP 


dx  hNAVULLyXrad 


(2.1-13) 


UL 


where  y  is  the  height  of  the  mirrors  or  active  medium,  whichever  is 
smaller. 

The  retention  of  the  spontaneous  emission  term,  S  ,  in  the  radia¬ 
tive  transfer  equations  has  resulted  in  an  explicit  expression  for  the 
intensity  in  terms  of  the  fluid  variables.  With  this  result,  the  Xj.aj 
terms  are  simply  nonlinear  terms  on  the  right  hand  sides  of  the  species 
equations.  Thresholds  and  cut  offs  are  obtained  as  part  of  the  solution 
of  the  differential  equations  for  the  flow  variables. 

The  lengthy  integration  times  associated  with  rotational  nonequili¬ 
brium  models  are  a  direct  consequence  of  the  number  of  dependent  variables 
and  the  stiffness  of  the  differential  equations.  While  the  number  of  de¬ 
pendent  variables  is  determined  by  the  number  of  levels  populated  by  the 
pumping  reaction,  the  stiffness  of  the  differential  equations  is 


i 
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determined  by  the  expression  for  I.  The  nonlinear  cascade  process  which 

leads  to  saturation  of  a  transition  is  primarily  controlled  by  the 

1-r  r  exp(2aL  )  term  in  the  denominator  of  Eq.  (2.1-11).  At  saturation, 

the  denominator  goes  to  zero  to  8  or  9  significant  figures;  the  resulting 

gain  curves  vs.  x  have  sharp  corners  at  saturation. 

A  reduction  factor  of  about  two  in  integration  time  can  be  realized 

by  using  the  approximate  expression  for  I,  which  was  derived  for  the  case 

aL  «  1 ,  Ref.  1 , 
e  ’ 


I  = 


4Kr  r .  aL 
_ o  L  e 

1-r  r.-2r  r.aL 
o  L  o  L  e 


(2.1-14) 


The  reduction  in  run  time  occurs  because  the  denominator  in  (2.1-14)  goes 
to  zero  linearly  rather  than  exponentially.  The  error  in  total  power  with 
this  approximation  is  less  than  0.6%.  The  only  problem  with  this  approxi¬ 
mation  is  that  it  does  not  give  the  correct  a  „  if  aL  is  not  much  less 

sat  e 

than  1,  which  is  the  case  for  many  high  power  chemical  lasers.  This  dif¬ 
ficulty  can  be  overcome  by  defining  an  effective  rQr^  which,  when  used  in 
the  approximate  expression  for  I,  gives  the  correct  a  .  This  is  ac- 

Sci  t 

complished  as  follows.  The  correct  asat  is  given  by 


sat 


-=r —  Jin  r  r 
2L  o  L 

e 


(2.1-15) 


From  the  linearized  equation  (2.1-14)  for  I,  the  approximate  a  is 

sat . 

A 


1-r  r. 
o  L 

lsat .  2r  r,  L 
A  o  L  e 


(2.1-16) 


By  setting  a  equal  to  a  ,  the  required  effective  r  r.  can  be  solved 
sat .  sat  o  l 

A  e 

for  from  Eq.  (2.1-16) 


o  I/eff  1+2L 

e  sat 


(2.1-17) 
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When  this  expression  for  (r0rj^eff  *s  uset*  *n  the  aPProximate  equation 

for  I,  Eq.  (2.1-14),  the  correct  a  will  be  obtained  from  the  approxi- 

sat 

e 

mate  I.  Thus,  the  present  model  employs  the  approximate  equation  for  I  with 
the  (r  r^)  ^  ^fitted  by  Eq.  (2.1-17). 

The  stiffness  of  the  differential  equations  can  be  controlled  by  the 
factor  K  which  is  the  frequency  average  of  the  spontaneous  emission  term. 
Increasing  K  decreases  the  stiffness  of  the  system  while  decreasing  K  in¬ 
creases  the  stiffness  of  the  system.  When  the  corners  of  the  gain  vs.  x 
curves  begin  to  round  off,  K  is  too  large.  The  factor  K  is  chosen  to  ensure 
that  the  corners  of  the  gain  vs.  a  curves  are  square.  Typically, 

K  =  ( 1 0_ 2  to  10)o'T‘*/tt  is  used. 

2.2  KINETICS  MODEL 

The  original  rotational  nonequilibrium  model  *  was  a  two  level,  pre¬ 
mixed  model  in  which  the  pumping  reaction  did  not  conserve  mass  and  deacti¬ 
vation  was  included  by  collision  with  a  generic  species.  To  obtain  quanti¬ 
tative  predictions,  the  correct  pumping  and  collisional  deactivation 
reactions  involving  the  primary  deactivators  HP,  F,  H  and  H2  have  been  in¬ 
cluded.  The  model  as  it  has  developed  to  date  comprises  the  following 


reaction  set: 

Pumping  Reactions 

F  +  H2  -  HF  (1 , J)  +  H  (2.2-1) 

F  +  112  ->  HF (2 ,J ' )  +  H  (2.2-2) 

F  +  H2  -  HF (3)  +  H  (2.2-3) 

F?  +  H  -*•  HF  (3)  +  F  (2.2-4) 


Collisional  Deactivation  Reactions 

HF(2,J')  +  M  Z  HF(1,J)  +  M  (2.2-5) 

HF(l.J)  +  M  t  HF COD  +  M  (2.2-6) 

HF (3)  +  M  t  HF(2,J')  +  M  (2.2-7) 

where  M  =  HF,  F,  H  and  H 

HF  (2 , J)  +  M  t  HF  (0)  +  M  (2.2-8) 

HF (3)  +  M  t  HF(1 ,J)  +  M  (2.2-9) 

HF (3)  *H?  HF (0)  +  M  (2.2-10) 

where  M  =  HF ,  F  and  H. 

Rotational  Relaxation  Reactions 

HF(2,J’)  +  M  t  HF(2,J)  +  M  (2.2-11) 

HF  (1 , J)  +  M  t  HF(1,J")  +  M  (2.2-12) 


where  M  =  HF,  F,  H,  F^,  He,  Ar. 

The  species  denoted  by  Ar  is  included  to  take  account  of  any  other  combus¬ 
tion  or  dissociation  products  that  may  be  present  in  the  mixture  and  which 
would  contribute  to  the  rotational  relaxation  but  not  to  the  collisional 
deactivation  of  the  lasing  species. 

While  the  above  model  is  more  extensive  than  originally  anticipated, 
it  is  quite  concise  compared  to  the  kinetic  models  employed  in  codes  such 
as  Blaze  11^  or  Lamp^.  To  duplicate  the  species  profiles  calculated  by 
Blaze  II,  it  was  necessary  to  add  the  reactions  (2.2-3),  (2.2-4),  (2.2-6), 
and  (2.2-7)  to  the  model.  The  addition  of  the  multiquantum  deactivation 
reactions  (2.2-8)  to  (2.2-10)  has  a  significant  effect  on  the  power  and 
the  length  of  the  lasing  zone,  see  Chapter  III. 
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The  rate  constants  are  input  in  the  usual  way  as 


n-N 


k  =  AT  exp(-B/RT) 


(2.2-13) 


where  A,  N  and  B  are  input  constants.  The  values  used  are  those  agreed 
upon  by  AFWL  and  Bell  Aerospace  Textron  in  June,  1977.  For  the  J-dependent 
reactions,  the  rate  constants  are  now  J-dependent  and  are  obtained  as 
follows.  For  the  pumping  reactions. 


k  =  kv  fv(J) 
pvJ  p  v  ’ 


(2.2-14) 


v 

where  k  are  the  global  forward  reaction  rates  which  have  been  determined 
P 

experimentally  as  a  function  of  temperature  and  which  are  used  in  the  ro¬ 
tational  equilibrium  models;  fV(J)  is  the  normalized  distribution  over  the 
rotational  states  of  each  v  level  (Z  fV(J)  =  1)  which  is  obtained  from  ex¬ 
periment.  For  HF  and  DF ,  the  data  of  Polanyi  and  Woodall^  is  used  for 
v 

f  (J).  The  pumping  reactions  are  assumed  to  be  one  way  reactions. 

For  the  collisional  deactivation  reaction  (2.2-5),  the  rate  of  change 
of  HF(2,J')  and  HF(1,J)  are  expressed  in  terms  of  the  deactivation  rate 
constants  as 


Xch  1 , J  P  ^  "m  j jL0  "HF  (2  ” J ' 5  k2  ,J  ;  1  ,J '  GJ ' 
+  J'=J+1  nHF(2’J,)  k2,J ' ; 1  ,J 


-  n^n.J)  exp(-hcw21/kT) 


y  kM 

L  T  .  1  it 

I  -n  *  > 1 


J'=0 


M 


,  M 


+  ^  k2  J ' • 1  J  ^J  I 

J'=J+1 


(2.2-15) 
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*c  h  2  ,J 


,  =  -  P2I 

M 


r"  tM  Jl 

+  J=jui  2.JJ1.J’  GJ'J 


-  exP(-hcw21/kT]  ^  nHF(1,J)k2,J' ;1,J  GJ 


l  n  (l,J)k' 


J=J'+1 


i  -i! 


(2.2-16) 


where  is  the  mole  mass  ratio  of  species  i.  Since  M  does  not  change 
during  the  reaction,  there  is  no  xch  term  for  the  M  in  these  reactions. 


Gj  =  (gj,/gj)  exp[- (EJ(  -  Ej)/kT] 


(2.2-17) 


gJ(  =  2J'  +  1,  gj  =  2J  +  1 


(2.2-18) 


Ej  =  hcBeJ(J  +  1) 


(2.2-19) 


■  2C1T-f'*1)e,p[-C2hc„:,1/l<T] 


(2.2-20) 


x  exp[-C2(EJt 


E^/kT] 


Z  .0>f«asxr/V 


aM  +  aHF 


(2.2-21) 


(2.2-22) 


(2.2-23) 


y  =  mHMHFMM/(MHF  +  V 


(2.2-24) 


and  Cn  are  input  constants  determined  by  requiring  that  the  following 


relation  be  satisfied 


k2,i  -J  zr 


(2.2-25) 


where  the  global  rate  constant,  which  is  determined  experimentally,  is 
expressed  as 

l  =  AT'Nexp(-B/RT)  (2.2-26) 

and  F(J')  is  the  Boltzmann  rotational  distribution  function.  When  M  is 
the  lasing  species,  n^  is 

nM  =  l  I  n(1F(v,J).  (2.2-27) 

v  J 

For  reaction  (2.2-6),  the  Xc^  terms  are 

Xch  =  -  P2["hf(1,J)  -  exp(-hcw1Q/kT)F(J)nHF(0)J 
1  ,J 

*  E  (2.2-28) 

«ch  ■  '  «P(-hc»10/kr)nl|F(o(]  E  vf*  (2.2-29) 

O  i—  — 1  M 

where  the  detailed  rate  constants  are  given  by 


.  M  M 

kl,J;0  =  ki,a 


(2.2-30) 


M  M 

k0;  1  ,J  =  k0,lF^ 


(2.2-31) 


where  F(J)  is  the  Boltzmann  rotational  distribution  function 


F (J)  =  exp[-  E  /kT] 

^rot 


Qrot  =  T/a9r 


(2.2-32) 

(2.2-33) 


and  a  is  1  for  a  heteronuclear  molecule  and  2  for  a  homonuclear  molecule. 


The  relations  (2.2-30)  and  (2.2-31)  between  the  detailed  and  the  global 
rate  constants  result  from  the  requirement  that  when  Eq.  (2.2-28)  is 


is  obtained.  When 


summed  over  J,  the  rotational  equilibrium  form  of  y  , 

ch 

M  is  the  lasing  species,  is  given  by  Eq.  (2.2-27). 

For  reaction  (2.2-7)  the  xch  terms  are  given  by 

Xch0  j,  =  p2[nHF(3)F(Jt)  -  nHF(2,J,)exp(-hca)32/kT)| 

X  J  VL  (2 . 2-34) 

Xch3  =  'p2[>(3)  ‘  nHF(2)exp(-hc(a32/kT)]  E  r^k”^  (2.2-35) 


where  the  detailed  rate  constants  are  related  to  the  global  rate  constants 
by 


k3-2  J'  =  k3  2F(J,:) 


k2,J';3  =  k2,3  ~  exp(-hcw32/kT)  ^ 


where  F(J')  is  the  Boltzmann  rotational  distribution  function. 

The  xc^  terms  for  the  multiquantum  deactivation  reactions 
(2.2-10)  are 


(2.2-36) 

(2.2-37) 


(2.2-8)  to 


Xch2  j ,  ~P  nMk2oj^HF^2,J'-)  '  nHF^0^exPC'haJ20/,kT^} 
r 

1  M 

Xch3  “p  ^k30^nHF^3^  ~  nHF  (^)exp  C"hcu)3Q/kT)  ] 

+  k31[nHF(3)  '  exP(-hcw31/kT)nHF(l)]  } 

XchQ  p  ^  |k30^nHF^3^  "  exP("hcw3o^kT^nHF 

+  k2ofnHF(2)  '  nHF(0)exP("hca)20/kT^ 

» 

Xch  =  p2nMk31[F(J)nHF(3)  '  exP('hcco31/kT)nHF(l,J)] 
1 ,  J 


(2.2-38) 


(2 . 2-39a) 


(2. 2-39b) 

(2.2-40) 


where  is  given  by  Eq.  (2.2-27)  when  M  =  HF  and  the  detailed  rate  constants 
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are  related  to  the  global  rate  constants  by 


M  M 

‘2 , J ' ;0  ~  k20 

(2 . 2-4  la) 

M  M 

0;2  ,J'  =F<J’>k20 

exp(-hcw20/kT) 

(2. 2-41b) 

'3;1  ,J  '  k31F(J) 

(2 . 2-42a) 

M  _  kM  _  kM 

1  ,J;3  13  '  31 

exp(-hcw31/kT) 

(2. 2-42b) 

Since  w20/T  ’  W31  /T  and  ^jq/T  are  greater  than  10  for  the  temperature  range  of 
interest  for  most  chemical  lasers,  the  backward  reaction  for  the  multi¬ 
quantum  deactivations  can  be  neglected.  With  this  approximation,  the  x„ 
terms  become 


xh 


(ch_,  j,  ^  "P  nHFt2,J,)  ^nMk20 


y  %  -~2 

Ch3 


p  n^t.!)  ^^30  +  k31^ 


Xcho  %  P2["hF^  >k30  +  "hF<2>  5"Mk20] 


XchUJ  -  P2nHF(^F^  >k31 


(2.2-43) 

(2.2-44a) 

(2. 2-44b) 

(2.2-45) 


where  F(J)  is  the  Boltzmann  rotational  distribution  function. 

For  the  rotational  relaxation  reactions  (2.2-11)  and  (2.2-12),  the 
Xcji  are  expressed  in  terms  of  the  deactivation  rate  constants  as 

fj-1 


Xchr  "p  Zr  (nHF  +^1nMEMr^ 


nnp  (i  »J) 


l  C 


J"=0 


M 


+  I 

J"=J+1 

JM 


X  pi»  <>  nHF^'j,,5 


l  Pj  nHF  (1  »J") 


J"=J+1 


(2.2-46) 
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X  chr  =  'P  Zr  (nHF  + 


HhpCZ.J') 


J'-l 


I  P 


J=0 


M 


i  p;J.  gj< 

j=j*+i  j  j 


j’-i 


J0  pj'  Cj'  "hf<2"» 


M 


J.J..1  PJ'  nHF(2,J) 


(2.2-47) 


where  nuc  =  Z  E  iwCv.J),  is  given  by  Eq.  (2.2-17),  the  rotational 

Hr  y  j  Hr  L 

transition  probability11  is  given  by 


PL  =  Cj  exp[-C2(Eu  -  EL)/kT] 


(2.2-48) 


Hp 

is  given  by  (2.2-22)  with  M  =  Hr .  Cj  and  C7  are  input  constants  de- 


.12,13 


termined  from  experiment*"’*''.  is  the  rotational  relaxation  efficiency 
of  species  M  compared  to  HF . 


2.3  FLUID  MECHANICS  MODEL 

To  keep  the  number  of  dependent  variables  to  a  minimum,  P(x),  T(x) , 
u(x) ,  the  mass  flow  rate  remaining  in  the  primary,  m^(x) ,  the  mass  flow 
rate  remaining  in  the  secondary,  ms(x),  and  the  ratio  of  the  thickness  of 
the  mixed  stream  to  the  geometric  thickness  of  the  flow,  Le/Lg,  are  input 
as  functions  of  x  from  a  detailed  rotational  equilibrium  model  such  as 
Blaze  II.  These  functions  are  input  as  polynomials  in  x  where  the  coef¬ 
ficients  are  determined  by  a  least  squares  fit  to  the  Blaze  II  profiles. 
Experience  has  shown  that  the  following  polynomials  give  the  best  fits  to 
these  variables: 


P  = 


C  Qx9  + 

p9 


C  «x' 

p8 


+  C  ,x’ 
P? 


+  C  ,x6 

p6 


+  C  x- 
P5 


♦  C  .x" 

p4 


C  ,x; 

p3 


C 

p2 


CplX 


+  C 


pO 


(2.3-1) 
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1  -  CT9x9  +  (-T8x8  +  CT7X?  +  cx6x6  +  CT5X$  +  CT4X^ 


*  CT3X  3  +  (T2x2  +  SlX  + 


(2.3-2) 


u  =  C  „x9+C  x8  +  C  7x7  +  C  x6  +  C  x5  +  C  .x4 
u9  u8  u/  u6  u5  u4 


♦  C  C  ^x2  ♦  C  .x  +  C 

uj  u2  ul  uO 


(2.3-3) 


m  =G1x+G_x1/2+mn 
P  Pi  P2  pO 


(2.3-4) 


m  =G,x  +  G7x1/2+iii. 
s  si  s2  sO 


Le/L«  =  V 


*  V,/2  *  V 


(2.3-5) 


(2.3-6) 


The  coefficients  of  these  polynomials  are  input  data.  For  the  ninth 

order  polynomials,  the  coefficient  must  be  input  to  nine  significant  figures. 

A  schematic  of  a  typical  mixing  region  is  shown  in  Fig.  2.  From  Ref.  6, 
the  mixing  term  which  appears  in  the  species  equation  is  given  by 


,  dm  dm 

M  .  =  —  (n.  -  n.  )  -j-~ -  +  (n.  -  n.  )  — 

ni  m  L  1  1S0  dx  1  ipo  dx 


(2.3-7) 


where  m  is  the  flow  rate  of  the  mixed  stream,  n^  is  the  mole  mass  ratio  of 

species  i,  n^so  is  the  value  of  the  mole  mass  ratio  of  species  i  in  the 

secondary  stream,  is  the  value  of  the  mole  mass  ratio  of  species  i  in 

the  primary  stream.  For  the  products  of  reaction,  n.  =  n.  =0;  for 

iso  lpo 

reactants  from  the  primary  stream,  n^sQ  =  0;  for  reactants  from  the 
secondary  stream,  n..  =  0.  The  mass  flow  rate  of  the  mixed  stream  is 

given  by 


m=m  -m  +  m  -m 
po  p  so  s 


(2.3-8) 


where  m  is  the  initial  primary  flow  rate,  msQ  is  the  initial  secondary 
flow  rate. 


Figure  2.  Schematic  of  the  mixing  region. 


20 


I  he  resulting  species  equations  are: 


dn  (1,J)  ; 

^ —  =  u  k;fi(j)nFnH. 


+  Z 


M  |Vc 


I  nHF(2’J')P2,J;l,J'GJ' 


J’=0 


+  j.=J+inHl:(2’J')P2’J,;1»J'exp('hcw2i/kT) 


+  y  pM  1 


nHF  ^ 


ZrF(nHF+^nMEMr) 


nHFC1’J) 


JV1  J  rM  J"  J" 
I  pj>>+  1  pj  G 

J"=0  J"=J+1  J  J 


A  PJ"GJ"nHF(1’J,,)  -  Jlt£J+1PJ  nHFC1’J") 


-  [nHF(l,J)-exp(-hcu,10/kT)F(J)nHFCO)]  E  n^ 


+  £  "Mk3!> 


+  ^  txCi.J.-i)  +  xU.J.D]  +  i 


dm  dm 

— EL  +  __s 
dx  dx 


nHF(1’J) 


(2.3-9) 


II 


dn  (2,J*)  ’ 

_ Mi _ £  V2f2 

dx  -  u  S,V 


+  ^  Pl  J-l  J'GJ'  'exP(-hcw  /kT)  £  nHr(l,J)p”  ,  ' 

J=J’+1  ’  ,I,J  J  21  a^Q  HF  J 


I  nHF(1’J)P2  J-l  J' 
J=J’+1  111 


:Cnl]F-E  n  E  )L;C2)J')fJi  f  P^, 

M  1  (j=0  J=Jr  +  l 


T 1  - 1 

^  r,  1  r  f  T  t 

V  r\°' 


PJ  GJ  nHF  ' J’'  “  J=j.  +  lPd’  nHF<'2,J'lJ  "  nHF(-2,J'^  ^nMk20 


+  [nHp(3)F(J)  -  nHF(2,J)exp(-hcu)32/kT)]  2j  r^k^ 


.  .  dm  dm 


(2.3-10) 


u  'ikPTnFnH„  “  kPHTnHnF^/ 


j  dm  dm 

+  m  nF  diT  +  (nF  "  nF0^  dx^ 


(2.3-11) 


c  |T>  9  C 


K  |  1  - c- 

“  PTnpnH2  *  •  (nH2  ■  nH2  >  sr  *  v  d/ 
L  20  2  J 


(2.3-1 2) 


fkPTnFnM2  ~  kPHTnHnF2J 


1 

~r~ 1 

dm 

dm 

1—5.  +  _ E 

m 

dx 

dx 

He 


poJ 


dm 
_£ 
dx  I 


' 

dm 

nA  “  nA 

Ar  Ar 

s 

dx  * 

so 

dm 

—2. 

dx 


(2.3-13) 


(2.3-14) 


2.3-15) 


2.3-16) 


.3-17) 


*  •  A  '  *  ■ 
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dnHF(3)  p 


dx  u  ]  nFnH  kp  "  tn||F(3)  *  nHF(2Jexp(-hcci»32/kT)]ZnMk32. 


nH*V  _,k PI l  '  nllF(-3^nM^k30  +  k31^ 


j  rdm  dm^- 

m  Ldx  +  dx  -  nMF^3^ 


(2.3-181 


where 


k 


PT 


k1  + 
P 


k3 

P 


(2.3-19) 


•'pHT  =  kpn  and  kpn  is  the  sum  of  kpn  over  v  =  3  to  6.  With  JM  =  20, 
there  are  50  differential  equations  for  the  mole  mass  ratios  of  the  various 
species  in  the  mixture.  The  initial  values  of  the  mole  mass  ratios  of 
the  species  are  calculated  from  the  initial  flow  rates  of  the  species 
in  the  primary  and  secondary  streams. 


m. 

n.  =  -22° 

lpo  m  W. 

po  1 


(2.3-20) 


n 


iso 


m  W . 
so  1 


(2.3-21) 


where  m  and  msQ  are  the  initial  primary  and  secondary  mass  flow  rates 

(m  =  E  m.  q  and  mso  =  Z  ) ,  ant*  w;  3S  t*ie  molecular  weight  of 
l  i  1  i 

species  i. 

In  some  cases,  F  atom  recombination  at  the  wall  is  important.  This 
effect  if  included  by  varying  the  F  atom  concentration  in  the  primary  stream 


/ 
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from  zero  at  the  wall  to  its  value  in  the  core  of  the  primary  flow.  This 


results  in  a  variation  of  n,.  ,  the  value  of  n.,  at  the  edge  of  the  mixing 

po 

layer  which  appears  in  the  mixing  terms,  with  x  as  shown  in  Fig.  3.  The 

boundary  layer  profile  is  used  to  provide  the  variation  of  n,.  from  the  wall 

r 

to  the  core.  The  mole  mass  ratio  is  related  to  the  species  mass  fraction  by 

ni  (y)  =  fj/W.  (2.3-22) 

o 

where  is  the  mass  fraction  of  species  i  and  ni  (y)  is  the  initial  value 
of  the  mole  mass  ratio  of  species  i  that  should  be  used  in  the  mixing  terms. 
The  value  of  y  depends  upon  the  value  of  x  and  the  thickness  of  the  mixing 
layer  y  =  6^  (x)  at  that  value  of  x.  Thus, 

M*'1 


nF  =  nF 


+  tnF  -  nF  ] 

wall  core  wall 


BL 


(2.3-23) 


:  =  ^ 

2  2wall 

o 


2c  ore 


2wall 


[Z.  ~1  3/2 


BL 


(2.3-24) 


Since  the  number  of  F  atoms  must  be  conserved,  np  is  related  to 

2wall 


nP  by 

wall 


1  1 

nF  2  nF  +  nF  ?  nF 

2wall  core  2core  wall 


(2.3-25) 


where  the  —  occurs  because  Wj;9  =  2Wp.  The  boundary  layer  profile  is  input 
from  Blaze  II  as 


<5,  00 


=  ax3/4  +  bx1/2 


(2.3-26) 


where  6^  =  BLT*bj  and  BLT  is  the  fraction  of  the  primary  nozzle  that  is 


I 


■ _ H — - 


.....  ,  Aw,,,.. 


stream  due  to  wall  recombination.  The  intersection 
of  the  F  atom  profile  with  the  edge  of  the  mixing 
region  provides  the  value  of  n_  (x)  which  is  used 

...  .  hpn 

in  the  mixing  term  in  the  species  equation. 
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filled  by  the  boundary  layer,  bj  is  the  half  height  of  the  primary  nozzle 
exit  plane,  and  a  and  b  are  curve  fit  constants  determined  by  fitting 
(2.3-26)  to  the  Blaze  II  data. 

The  final  element  of  the  model  is  the  inclusion  of  the  effect  of 
lasing  on  the  1  to  0  band.  This  was  necessary  to  obtain  quantitative 
agreement  with  the  Blaze  II  results.  Since  the  present  model  is  a  two  level 
model,  the  effect  of  lasing  on  the  1  to  0  band  was  included  as  a  rate  con¬ 
stant  which  was  determined  from  the  Blaze  II  value  of  xj.^00  by 

setting 


X 


1  o 
rad 


krad  P  nHFnHF('1-) 


(2.3-27) 


which  simulates  the  depopulation  of  HF(1)  by  lasing  as  a  collisional  de¬ 
activation  by  HF .  From  Eq.  (2.3-27),  (x)  is  obtained  and  input  as  a 
function  of  x  in  the  program.  The  rate  constant  for  the  deactivation  of 
HF(1)  by  HF  is  then  written  as 


k 


10 


k10 

rad 


(2.3-28) 


and  the  effect  of  lasing  on  the  depopulation  of  HF(1)  is  taken  into  account 
in  the  collisional  deactivation  reaction. 


Since  the  present  model  is  a  two  level  model  as  far  as  the  calculation 
of  laser  power  is  concerned,  the  gain  is  given  by 


a(v,J,m)  =  PhNAwUL<J>CwUL)BL(J[(gL/gu)nu  -  nj  (2.3-29) 


where  for  a  doppler  profile 


V*D  =  C 


!k_ 

2ttR 


1/2 


p'  1  /z 


(2.3-30) 
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WL  is  the  molecular  weight  of  the  lasing  species. 


(2.3-31) 


gL  =  2J  +  1’  gU  =  2J  +  1  +  2m  (2.3-32) 

and  m  =  -1  for  the  P  branch  and  m  =  +1  for  the  R  branch.  For  a  Voight 
profi le , 

^Voight  =  VC°’y)  (2.3-33) 

where  K(o,y)  is  given  in  Ref.  6.  The  power  in  each  transition  is  given  by 

hN  cw..  i-x 

P(x,v,J,m)  =  — ^ -  A  TX(v,J,m)dx  (2.3-34) 

JM  Jo 

where  AmT  =  m/pu,  m  is  the  total  flow  rate  in  the  mixed  flow  and 

FJm  =  a)21/u)UL  =  +  [m(2J+l+m)erk/hcu21]}'1  (2.3-35) 

The  density  p  is  given  by 

p  =  PW/RT  (2.3-36) 

where  the  molecular  weight  of  the  mixture  is  given  by 


W  =  (In.)-1 
i  1 


(2.3-37) 


and  R  is  the  universal  gas  constant  per  mole. 

In  the  next  chapter,  the  present  model,  which  is  denoted  by  the  name 
MNORO,  is  compared  to  the  Blaze  II  results  for  the  arc  driven  laser  as  well 
as  for  the  combustion  driven  CL  XI  case.  From  these  comparisons,  the 
essential  elements  of  a  quantitative  model  can  be  identified. 


....  . .. 


.K 
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III.  COMPARISON  WITH  THE  BLAZE  II  MODEL 
In  attempting  to  develop  a  simplified  cw  chemical  laser  model  that  will 
give  quantitative  predictions,  it  is  necessary  to  carry  out  as  detailed  com¬ 
parisons  with  other  models  or  experimental  data  as  possible.  In  the  pre¬ 
sent  study,  two  diverse  cases  were  selected  for  initial  baselining  of  the 

9 

model.  One  is  an  Aerospace  arc  driven  cw  HF  power  spectral  case  and  the 

14 

other  is  the  Bell  Aerospace  Textron  CL  XI  lasing  lines  case 

3.1  THE  AEROSPACE  ARC  DRIVEN  LASER 
The  Aerospace  nozzle  is  shown  in  Fig.  4.  The  case  chosen  for  base- 

9 

lining  was  the  power  spectral  run  ,  the  flow  conditions  for  which  are  given 
in  Table  I.  The  optical  cavity  was  a  Fabry-Perot  resonator  composed  of 
two  flats.  The  spectra  was  taken  by  hole  outcoupling  through  one  of  the 
flat  mirrors  at  1.905  cm  downstream  from  the  nozzle  exit  plane.  In  both 
the  Blaze  II  and  the  present  model,  the  SF&  was  assumed  completely  dis¬ 
sociated;  the  sulfur  and  0^  were  treated  as  inert  species  which  were 
represented  by  Argon  in  both  models.  The  experimental  power  spectral 
distribution  is  shown  in  Fig.  5.  Because  of  the  uncertainties  in  the 
modeling  of  the  spray  bar  mixing  and  in  the  percent  dissociation  of  SF,  in 

D 

the  arc,  no  real  attempt  to  match  the  experimental  data  was  made;  rather, 
the  elements  in  the  present  simplified  model  which  are  required  to  match 
the  Blaze  II  results  for  the  same  initial  conditions  are  identified  through 
the  comparison  of  the  two  models. 

The  Blaze  II  fluid  dynamic  profiles  together  with  the  resulting  poly¬ 
nomial  fits  used  in  MN0R0  are  shown  in  Figs.  6-10.  As  can  be  seen,  the 
polynomials  give  an  excellent  representation  of  the  fluid  dynamic  variables. 

The  Blaze  II  and  MNORO  H,  F  and  H^  profiles  are  compared  in  Fig.  11. 

As  long  as  the  production  of  HF(3)  by  the  pumping  reaction  is  included. 


MULTIPLE  N02ZLE 


Figure  4.  Schematic  of  the  Aerospace  laser  nozzle. 
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Species 

in  gm/s 

L,  cm 

SF6 

0.5 

17.78 

He 

2.06 

°2 

0.75 

T  °K 

P 

H2 

1.25 

160 

Table  I.  Aerospace  CIV  HF  Laser  Power  Spectral  Run  Conditions. 


/ 


P  torr 

1.4 

T  °K 
s 

80 


BLAZE  It 
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Aerospace  laser  velocity  profile. 
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Aerospace  laser  mass  flow  rate  remaining  in  the  secondary  stream  (m  )  profile 


36 


O 

H 


AXIAL  DISTANCE,  CM 

Figure  10.  Aerospace  laser  effective  length  of  the  active  media  divided  by  the  geometric  length  of 
the  flow  channel  in  the  lasing  direction  as  a  function  of  x. 


AL  SPECIES  PROFILES 
-  BLAZE  31 


excellent  agreement  with  the  Blaze  II  H,  F  and  Il2  profiles  is  obtained 
regardless  of  whether  lasing  on  3  2  is  simulated  by  setting  kp2  =  kp2+kp3 

or  whether  HF(3)  is  pumped  and  allowed  to  collisionally  deactivate.  This 
agreement  is  independent  of  whether  the  model  is  run  in  the  rotational 
equilibrium  or  the  rotational  nonequilibrium  mode. 

The  Blaze  II  and  MNORO  HF(v)  (v  =  0,  1,  2,  3)  profiles  are  shown  for 
the  rotational  equilibrium  case  in  Fig.  12  and  for  the  rotational  nonequi¬ 
librium  case  in  Fig.  13.  Since  Blaze  II  has  14%  of  the  total  power  in  the 
3  -*•  2  band,  we  see  that  the  best  agreement  with  the  HF(v)  profiles  occurs 
when  lasing  on  3  +  2  is  simulated  by  setting  kp2  =  kp2+kp3.  When  lasing  on 
3  -*■  2  is  not  simulated,  we  see  that  multiquantum  deactivation  has  a  neg¬ 
ligible  effect  on  all  the  HF(v)  profiles  except  for  HF(3)  which  shows  a 
small  effect.  When  lasing  on  3  +  2  is  not  included,  the  rotational  non¬ 
equilibrium  profiles  tend  to  deviate  more  from  the  Blaze  II  profiles  than 
the  equilibrium  ones  do.  The  fact  that  in  the  rotational  nonequilibrium 
case  the  lasing  deactivation  of  HF (2)  is  slower  than  in  the  equilibrium 
case  results  in  more  HF(2)  and  less  HF(1).  The  significant  effect  of  lasing 
on  3  +  2  on  the  HF(3)  profile  is  evident  in  Fig.  12  by  comparing  the  Blaze 
II,  which  includes  lasing  on  3  -*■  2,  and  the  MNORO  case  for  pumping  HF(3)  and 
subsequent  collisional  deactivation. 

The  power  spectral  distributions  are  compared  in  the  equilibrium,  Fig.  14, 
and  the  nonequilibrium.  Fig.  15,  cases.  As  can  be  seen,  lasing  on  3  2  and 

the  different  types  of  collisional  deactivation  have  only  a  small  effect  on 
the  2  ->•  1  power  spectral  distribution. 

From  Table  II,  it  is  seen  that  when  lasing  on  3  +  2  is  simulated  in 
MNORO,  the  power  agrees  very  well  with  the  Blaze  II  power.  Thus,  the  MNORO 
model  gives  very  good  agreement  with  the  Blaze  II  species  profiles,  power 
and  power  spectral  distribution. 
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MNORO  SIMULATE  3-2  LASING  BY  SETTING 


Figure  13.  Comparison  of  the  Blaze  II  and  the  rotational  nonequilibrium 
MNORO  HF(v)  profiles  for  the  Aerospace  laser  (AL) . 


Comparison  of  the  Blaze  II  and  MNORO  power  spectral  distributions  for  the  Aerospace  laser  CAL) 


15.  MNORO  rotational  nonequi librium  power  spectral  distributions  for  the  Aerospace  laser  (AL) . 


BLAZE  II 


485 


515 


1 


MNORO  ROTATIONAL  EQUILIBRIUM 

Simulate 

529 

3-2  Lasing 

Lasing  on  2-1,  1-0 

384 

Lasing  on  2-1,  1-0,  MQVT 

375 

MNORO  ROTATIONAL  NONEQUILIBRIUM 


Simulate 

500 

3-2  Lasing 

Lasing  on  2-1,  1-0 

362 

Lasing  on  2-1,  1-0,  MQVT 

347 

Table  II.  Comparison  of  BLAZE  II  AND  MNORO  2  +  1  P. 
for  the  Aerospace  Laser. 
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Several  other  effects  are  evident  in  Table  II.  First,  when  mixing  is 
present,  the  rotational  nonequilibrium  power  is  only  about  6%  less  than  the 
rotational  equilibrium  power  rather  than  the  20%  difference*  obtained  in  a 
premixed  case.  Second,  lasing  on  3  +  2  has  a  significant  effect  on  the  2  -+  1 
power.  When  14%  of  the  total  power  is  in  the  3  ■+  2  band,  the  2  -*•  1  power 
increases  38%  in  both  the  equilibrium  and  the  nonequilibrium  cases.  This 
suggests  that  it  might  be  worthwhile  to  promote  3  -*•  2  lasing  in  order  to  en¬ 
hance  the  2  -*■  1  power.  On  the  other  hand,  it  points  out  the  danger  of  ad¬ 
justing  parameters  in  a  model  to  produce  agreement  with  total  power  if  the 
model  predicts  lasing  on  3  -*•  2  when  it  is  not  observed  experimentally.  These 
results  suggest  that  there  may  be  a  problem  with  the  kinetic  model.  We  will 
see  this  same  effect  in  the  CL  XI  case,  but  even  more  pronounced  there,  since 
Blaze  II  predicts  lasing  on  the  upper  bands  as  well  as  3  +  2,  none  of  which 
is  observed  experimentally. 

3.2  THE  BELL  AEROSPACE  TEXTRON  CL  XI  EXPERIMENT 

Bell  Aerospace  Textron  has  performed  a  series  of  diagnostic  experiments 
14 

with  the  CL  XI  nozzle  .  The  CL  XI  (15,30)  HF  lasing  lines  experiment  was 
chosen  for  initial  baselining  of  the  present  model.  In  this  experiment,  the 
power  spectral  distribution  and  the  intensity  distribution  on  each  line 
were  measured.  The  intensity  distributions  were  measured  with  cylindrical 
mirrors  oriented  with  the  axis  of  the  cylinder  parallel  to  the  x  axis  (flow 
direction).  Thus,  a  Fabry-Perot  model  of  the  cavity  should  provide  a  good 
representation  of  this  experiment. 

The  primary  flow  contained  F,  He,  N^,  CF^  and  DF  and  the  secondary 
flow  contained  and  He.  The  N^,  CF^  and  DF  were  treated  as  inert  species 
which  were  represented  by  Argon  in  the  present  simplified  model.  To  pre¬ 
serve  numbers  of  particles,  the  total  molar  flow  rate  of  these  species  was 
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replaced  by  the  same  molar  flow  rate  of  Argon.  The  experimental  power 
spectral  and  intensity  distributions  are  shown  in  Figs.  16  and  17.  As  in 
the  previous  case,  the  emphasis  to  date  has  been  on  comparing  the  present 
MNORO  model  run  in  the  equilibrium  mode  with  the  Blaze  II  results  for  the 
CL  XI  (15,30)  case. 

Several  Blaze  II  runs  were  made  for  this  study.  The  first  one  allowed 
all  bands  up  to  8  -*■  7  to  lase.  The  second  case  allowed  only  3  2 ,  2  ■+■  1 

and  1  -►  0  to  lase.  The  last  case  allowed  only  2  -*■  1  and  1  -*•  0,  the  bands 
observed  experimentally,  to  lase. 

The  Blaze  II  calculations  took  the  effect  of  the  nozzle  trip  flow  on 
the  mixing  into  account  by  increasing  the  binary  diffusion  coefficient. 

This  multiplier  was  established  by  comparison  of  lasing  zone  length  and 
total  power  with  experimental  data  for  the  case  when  all  bands  up  to  8  -*■  7 
could  lase  if  conditions  would  allow  it1'’.  The  results  of  these  three  runs 
are  shown  in  Table  III.  The  effect  of  the  change  in  the* vibrational  bands 
which  were  allowed  to  lase  on  the  fluid  dynamic  and  species  profiles  is  shown 
in  Figs.  18-20.  From  these  figures,  it  is  seen  that  the  hot  reaction  bands 
have  a  negligible  effect  on  these  profiles.  Not  allowing  lasing  to  occur  on 
3  2  results  in  a  10%  increase  in  the  peak  pressure  and  temperature  and  a 

negligible  effect  on  the  other  profiles  except  for  the  HF (3)  profile.  Thus, 
one  set  of  fluid  profiles  was  used  throughout  this  study.  As  in  the  preceed- 
ing  case,  the  polynomials  give  an  excellent  representation  of  the  fluid  dyna¬ 
mic  variables. 

Figure  21  shows  the  effect  of  the  bands  which  are  allowed  to  lase  on  the 
2  ->  1  power  spectral  distribution.  The  hot  reaction  bands  have  a  negligible 
influence  on  the  2  -*■  1  power  spectral  distribution.  However, 


P,CJ)/PT  B,(J)/P 


Bands  Lasing  Lasing  Zone  Total 


Figure  18.  Blaze  II  CL  XI  (15,30)  temperature  and  pressure  profiles 
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Figure  21.  Blaze  II  CL  XI  (15,30)  power  spectral  distributions. 
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lasing  on  3  -►  2  has  a  significant  effect  not  only  on  the  total  power  in  the 
2  *►  1  band.  Table  III,  but  also  on  the  2  -►  1  power  spectral  distribution. 

The  MNORO  results  excluding  the  hot  reaction  are  compared  to  the  Blaze 
II  species  profiles  in  Figs.  22  and  23  for  the  cases  of  no  lasing  on  3  +  2 
and  lasing  on  3  -*■  2  simulated  by  setting  kp2  =  kp2  +  kp3.  The  corresponding 
power  spectral  distributions  and  2  -*■  1  powers  are  compared  in  Fig.  24.  From 
these  figures  we  see  that  even  though  the  P-^'s  agree  within  10%  and  the 
power  spectral  distributions  are  quite  similar,  there  are  considerable  dif¬ 
ferences  in  all  the  species  profiles  except  for  the  H  atom  profile. 

The  corresponding  results  when  the  hot  reaction  is  included  by  pumping 
only  to  HF(3)  and  when  multiquantum  VT  deactivation  by  I1F,  F  and  H  are  in¬ 
cluded  are  shown  in  Figs.  25-27.  Because  of  the  relatively  high  temperature 
(T,MAX  =  585°K)  and  the  F  atom  wall  recombination  which  results  in  a  signi¬ 
ficant  amount  of  F?,  the  hot  pumping  reaction  must  be  included  in  the  kinetic 
model  to  duplicate  the  Blaze  II  F,  H  and  F?  profiles.  Fig.  25.  As  might  be 
expected,  the  multiquantum  deactivation  has  a  negligible  effect  on  the  F,  H, 

F  and  H_  species  profiles.  From  Fig.  26,  it  is  seen  that  both  the  hot  re- 
action  and  multiquantum  deactivation  are  required  to  obtain  good  agreement 
with  the  Blaze  II  HF(v)  profiles.  The  agreement  is  not  as  good  for  HF(3) 
because  the  hot  reaction  is  pumping  completely  to  HF(3)  rather  than  to  the 

upper  levels  with  collisional  cascade  to  HF (3) .  When  hot  reaction  pumping 

6 

l  n  y 

to  HF(3)  and  HF  (4) ,  where  kp^  =  \  inclu<led>  there  is  no  change  in 

v=4 

the  HF(v)  species  profiles  for  v  =  0,  1,  2  and  only  a  slight  improvement  in 
HF(3).  To  obtain  better  agreement  with  HF(3)  would  require  including  at 
least  up  to  HF(5)  because,  based  on  the  magnitude  of  the  rate  constants  in¬ 
volved,  there  is  a  bottle  neck  at  HF(4)  which  would  result  in  decreased 
HF(3).  Since  the  effect  on  the  2  1  power  and  power  spectral  distribution 
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Figure  23.  Comparison  of  Blaze  II  and  MNORO  HF(v)  profiles  for  the  CL  XI  (15,30), 
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Figure  24.  Comparison  of  Blaze  II  and  MNORO  power  spectral  distributions  for  CL  XI  (15,30). 
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Figure  26.  Compvirison  of  Blaze  II  and  MNORO  HF(v)  profiles  for  the  CL  XI  (15,30) 
showing  the  effect  of  the  hot  reaction  and  MQVT. 
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would  be  small,  the  additional  complexity  is  not  warranted.  Thus,  the  final 

model  includes  the  effect  of  the  hot  reaction  by  allowing  it  to  pump  to 

6 

HF(3)  only,  i.e.,  The  e^^ect:  °f  the  hot  reaction  and  multi- 

v=3 

quantum  deactivation  on  the  rotational  equilibrium  power  spectral  distribu¬ 
tion  is  shown  in  Fig.  27.  Again  it  is  seen  that  both  of  these  reaction  types 
have  a  significant  effect  and  hence  must  be  included  in  the  model.  From  runs 
made  with  only  the  2  -*■  0  multiquantum  deactivation  and  with  2  -*•  0,  3  -*■  1  and 
3  ■+•  0,  the  most  important  of  these  reactions  is  the  2  -*■  0  deactivation  by  a 
factor  of  about  two  to  one. 

The  rotational  nonequilibrium  species  profiles  when  the  hot  reaction 
pumps  only  to  HF(3)  and  multiquantum  deactivations  are  included  are  compared 
to  the  Blaze  II  profiles  in  Fig.  28  and  29.  As  can  be  seen,  the  nonequili¬ 
brium  species  profiles  agree  as  well  with  the  Blaze  II  species  profiles  as 
the  equilibrium  ones  do.  However,  it  is  seen  from  Fig.  30  that  the  nonequi¬ 
librium  power  spectral  distribution  is  considerably  different  than  the 
equilibrium  power  spectral  distribution.  The  peak  of  the  nonequilibrium  PSD 
at  J  =  8  is  indicative  of  the  relatively  large  T  (T  =  585°K)  :  since  the 
experimental  data  peaks  at  J  =6,  Blaze  II  may  be  overpredicting  the  tem¬ 
perature. 

Figures  31-36  give  the  rotational  equilibrium  and  the  rotational  non¬ 
equilibrium  saturated  gains  and  intensity  distributions.  The  J-shifting  is 
evident  in  the  gain  and  intensity  distributions  for  the  equilibrium  case 
whereas  simultaneous  lasing  on  many  lines  occurs  in  the  nonequilibrium  case. 

From  the  Blaze  II  results.  Table  III,  it  is  seen  that  the  occurrence 
of  lasing  on  3  ->  2  results  in  a  21%  increase  in  the  power  on  2  +  1;  MNORO 
indicates  a  24%  increase  when  lasing  on  3  -*■  2  occurs.  Since  this  is  a  sig¬ 
nificant  effect,  as  found  in  the  previous  case,  it  may  be  worthwhile 
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CL  XI  SPECIES  PROFILES 


Figure  28.  Comparison  of  Blaze  II  and  rotational  nonequilibrium  MNORO 
II,  F,  H2  and  F2  profiles  for  CL  XI  (15,30). 
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XXX  MNORO  ROTATIONAL  NONEQUILIBRIUM 

HOT  REACTION  PUMPS  HF(3),MQVT 
LASING  ON  2-1, 1-0 


Figure  29.  Comparison  of  Blaze  II  and  rotational  nonequilibrium 
MNORO  HF (v)  profiles  for  CL  XI  (15,30). 


Comparison  of  rotational  equilibrium  (RE)  and  rotational  nonequilibrium  (RNE) 
MNORO  power  spectral  distributions. 
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promoting  3  -*  2  lasing  to  obtain  the  increased  performance  on  2  -»  1  and  1  ♦  0. 
As  in  the  previous  case,  the  equilibrium  and  the  nonequilibrium  powers  only 
differ  by  about  4%  rather  than  the  20%  difference1  observed  in  premixed  cases. 

From  these  comparisons,  it  is  seen  that  the  MNORO  model  includes  the 
essential  elements  required  to  give  good  agreement  with  the  very  detailed 
Blaze  II  model.  These  elements  are:  lasing  on  1  +  0  and  3  -*■  2  if  it  occurs, 
F-atom  wall  recombination  and  the  hot  reaction.  The  importance  of  the  multi¬ 
quantum  deactivation  depends  upon  the  temperature.  For  the  Aerospace  device, 
in  which  the  temperature  was  quite  low  (n/230°K)  ,  multiquantum  effects  were 
less  than  3%;  however,  for  the  CL  XI  device,  in  which  the  temperature  was 
quite  high  (>5850K) ,  multiquantum  deactivation  results  in  about  a  29%  de¬ 
crease  in  power  and  has  a  significant  effect  on  the  length  of  the  lasing 
zone.  Table  IV. 

This  section  will  be  concluded  with  a  few  remarks  on  run  time.  The  ro¬ 
tational  nonequilibrum  model  employs  sparse  matrix  techniques  in  the  evalua¬ 
tion  and  inversion  of  the  Jacobian  matrix  used  by  the  implicit  integration 
scheme.  The  use  of  sparse  matrix  techniques  reduces  run  time  by  factor  of 
three.  A  typical  case  with  the  original  rotational  nonequilibrium  model, 

Noro  I,  requires  30  sec  (all  times  are  on  the  CYBER  175).  The  MNORO  model 
varied  from  108  sec  for  the  Aerospace  case  to  198  sec  for  the  CL  XI  case. 

At  this  time,  the  reason  for  the  factor  of  2  increase  from  one  case  to  the 
next  is  not  clear.  It  may  be  associated  with  the  fact  that  the  Aerospace 
case  is  a  slower  mixing  device  and  thus  the  physics  occurs  over  longer  dis¬ 
tances  permitting  larger  integration  steps.  This  is  suggested  by  the  fact 
that  the  Aerospace  case  took  1601  steps  wherease  the  CL  XI  case  took  3123 
steps  to  go  10  cm.  It  should  be  kept  in  mind,  however,  that  the  effect  of 
the  temperature  profiles  on  run  time  is  not  clear;  the  Aerospace  case  had  an 
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Model 


Blaze  II 
MNORO 

NO  Hot  Reaction 
NO  MQVT 

Hot  Reaction 
NO  MQVT 

Hot  Reaction 
MQVT 

ROTATIONAL 
Nonequilibrium 
Hot  Reaction 
MQVT 


^21  *Cut  off 

1.0  1.0 


0.89 

1.27 

0.98 


1.87 

1.57 

1.15 


0.93 


1.25 


Table  IV.  The  Effect  of  the  Hot  Reaction  and  Multiquantum  VT  Deactivation 
on  the  2  1  Power  and  the  Length  of  the  Lasing  Zone  for 

CL  XI  (15,  30). 


72 


almost  constant  T  whereas  the  CL  XI  profile  is  similar  to  an  inverted  N,  and 
since  the  kinetic  rates  are  exponential  in  T,  this  difference  may  be  signifi¬ 
cant.  The  basic  increase  in  run  time  from  30  sec  for  the  original  model  to 
108  sec  for  MNORO  is  a  consequence  of  the  increase  in  the  number  of  Xc^em 
terms  that  must  be  evaluated.  The  original  model  had  only  one  deactivating 
species  whereas  MNORO  has  four.  The  addition  of  dependent  variables  that  are 
not  lasing  species  affects  run  time  only  through  the  number  of  additional 


^chem  terms  to  be  evaluated.  The  run  time  depends  quite  strongly  on  the 
stiffness  of  the  differential  equations,  which  is  determined  by  the 
terms.  When  a  transition  lases,  the  number  of  places  to  which  a  agrees  with 


-  — )  £nr  r. 

vLe'  o  L 


is  critical.  This  can  be  controlled  by  the  quantity  K 


(see  Section  2.1).  Generally,  using  K  =  10  o'T4/tt  results  in  a  agreeing 


with  a  to  5  or  6  places  which  requires  the  integration  accuracy  parameter 
e  =  10  6.  This  combination  of  K  and  e  appears  to  minimize  run  time  with  no 


loss  of  accuracy. 


Finally,  it  should  be  pointed  out  that  the  MNORO  model  cannot  be  run 
in  the  equilibrium  mode  by  increasing  the  rotational  relaxation  rate  constant 
by  a  factor  of  105  when  the  integrator  employs  sparse  matrix  techniques.  To 
operate  MNORO  in  the  quilibrium  mode  requires  the  complete  Jacobian  which 
generally  results  in  a  factor  of  5  increase  in  run  time. 

From  the  preceeding  results,  it  is  seen  that  the  MNORO  model  developed 
in  this  study  is  an  efficient  rotational  nonequilibrium  model  which  can  be 
expected  to  provide  quantitative  predictions  of  the  power  and  power  spectral 
distribution  of  cw  chemical  lasers. 


I 


0 


73 


IV.  COUPLED  PHYSICAL  OPTICS  .AND  KINETICS 
The  second  phase  of  the  study  involved  the  coupling  of  the  AFWL  strip 
resonator  physical  optics  code  supplied  by  Captain  Salvi,  denoted  by  AFOPT, 
with  the  rotational  nonequilibrium  kinetics  model.  This  code  has  been 
coupled  to  the  original  rotational  nonequilibrium  model1,  denoted  N0R0-1, 
and  to  the  model  developed  during  this  program,  denoted  MNORO. 

In  coupling  the  kinetics  to  the  physical  optics,  the  gain  medium  is 
treated  as  a  thin  sheet  located  between  the  mirrors.  The  AFOPT  code  has 
the  capability  of  allowing  an  arbitrary  placement  of  the  gain  sheet  by 
changing  the  sequence  in  which  various  subroutines  are  called.  To  check 
that  the  coupling  to  the  kinetics  was  correct,  several  previous  cases'’, 
which  used  the  Bell  Aerospace  Textron  strip  resonator  code,  denoted  BATOPT, 
were  repeated  with  the  AFOPT  +  NORO-I  code.  Since  these  previous  calcula¬ 
tions  placed  the  gain  sheet  adjacent  to  the  large  mirror,  this  placement 
of  the  gain  sheet  was  used  here.  The  geometry  of  the  unstable  confocal 
resonator  is  shown  in  Fig.  37.  The  intensity  sent  to  the  kinetics  must  be 
the  average  of  the  left  and  right  running  waves, 

JAV  =  ^ISt1+exp(°tL8)]  +  I3(l  +  exp(aLg)]  }/2  (4-1) 

where  I<.  and  I^  are  the  intensities  at  positions  S  and  3  respectively.  The 
solution  is  an  iterative  one.  The  calculation  starts  by  passing  a  plane 
wave  of  roughly  the  expected  intensity  through  the  empty  cavity.  The  re¬ 
sulting  intensity  distribution  is  sent  to  the  kinetics.  The  species  equa¬ 
tions  are  solved  with  the  given  intensity  distribution,  and  the  gains  on 
each  line  are  calculated  and  sent  to  the  optics.  The  stored  wave  from  the 
preceeding  optics  calculation  is  then  propagated  through  the  cavity  and  in¬ 
cremented  by  the  gain  distribution  from  the  kinetics.  This  process  is  re¬ 
peated  until  the  difference  between  successive  iterates  is  less  than  some 
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prescribed  amount.  Generally,  the  error  continues  to  decrease  as  the  num¬ 
ber  of  iterations  increases. 

To  check  the  convergence  of  the  iteration  procedure,  it  has  been  found 
useful  to  compare  the  intensity  distributions  of  successive  iterates  and 
to  compare  the  power  extracted  from  the  medium  as  calculated  in  the  kinetics 
and  the  optics  parts  of  the  code.  The  power  calculated  in  the  kinetics 
part,  pchem»  is  the  total  power  extracted,  both  outcoupled  and  the  dif- 
fractively  lost  power;  thus,  pchom  must  be  compared  to  pLQSS  TOtai  »  the 
total  power  lost  from  the  resonator  in  the  optics  part  of  the  code.  Re¬ 
ferring  to  Fig.  37^ 


P  =P-P+P-P 

LOSS  TOTAL  1  3  *4  5 


(4-2) 


where  P^  is  the  power  at  station  i.  In  addition  to  the  calculation  of  P^, 

P1  ’  P2’  P2a’  P3’  P4  an^  P5’  a  Printer  piQt  of  the  intensity  distribution 
at  the  outcoupling  mirror,  location  2,  for  each  line  and  for  the  sum  over 
all  lines  was  added  to  AFOPT.  This  printer  plot  is  quite  useful  for  moni¬ 
toring  the  changes  between  successive  iterates. 

The  first  case  run  with  the  AFOPT  +  NORO-I  code  was  the  50%  out- 
coupled,  confocal  unstable  resonator  for  the  single  line  J  =  8^.  This 
case  was  selected  because  it  converged  in  15  iterations  with  no  oscilla¬ 
tions.  The  first  three  iterations  were  run  several  times  with  various  com¬ 
binations  of  the  AFOPT  control  parameters.  Since  quite  different  intensity 
distributions  were  obtained  depending  upon  the  choice  of  the  parameter  grid- 
rad,  an  empty  resonator  case  was  run  to  convergence  to  study  the  role  of  the 
guard  band  on  the  convergence  and  accuracy.  The  empty  cavity  case  was  used 
because  it  minimized  CPU  time. 

The  empty  cavity  solution  for  a  confocal  unstable  resonator  with  50% 
geometric  outcoupling  was  obtained  with  the  AFOPT  and  BATOPT  codes.  At 
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convergence,  the  two  solutions  agreed.  However,  particularly  at  early 
iterations,  there  were  considerable  differences  between  the  two  codes. 

These  differences  were  perplexing  since  the  two  codes  purported  to  follow 
the  same  physical  process  using  identical  initial  conditions. 

It  was  theorized  that  the  observed  discrepancies  might  be  associated 
with  the  spatial  frequency  content  of  the  field  distribution.  During  the 
first  few  round  trips  through  the  cavity,  the  spatial  frequency  content  or 
nonuniformity  of  the  field  distribution  would  be  higher  than  in  later  tran¬ 
sits.  As  convergence  was  approached,  the  field  would  become  increasingly 
more  uniform  as  the  lossy  high  frequency  components  diffractively  spread 
out  of  the  resonator.  Thus,  the  observed  difference  might  be  ascribed  to 
a  variation  in  the  effectiveness  of  the  high  frequency  diffractive  pro¬ 
pagators  employed  in  the  two  codes.  This  would  be  a  critical  factor  af¬ 
fecting  the  outcome  of  the  overall  study  since  the  addition  of  gain, which 
only  parially  fills  the  resonator,  stimulates  and  maintains  the  higher 
spatial  frequency  components.  This  rationale  was  the  basis  for  the  con¬ 
tinued  intensive  investigation  of  the  early  iterate  anomalies. 

As  background,  the  major  similarities  and  differences  between  the  two 
codes  will  be  reviewed.  Both  codes  begin  with  an  arbitrary  input  distri¬ 
bution  located  at  some  point  in  the  resonator  and  diffractively  propagate 
it  through  the  resonator  until  a  stable  irradiance  and  phase  distribution 
are  obtained.  The  version  of  the  Bell  Aerospace  Textron  ID  resonator  code 
being  used  in  this  case  always  expresses  the  complex  field  on  a  plane  per¬ 
pendicular  to  the  optical  axis.  When  mirrors  are  encountered,  the  phase 
of  the  field  on  the  plane  is  adjusted  to  approximate  the  effects  of  the 
curvature  of  the  mirror.  The  amplitude  of  the  field  at  a  mirror  is  ad¬ 
justed  to  correct  for  the  mirror  absorption  over  the  clear  aperture  and 
the  edges  are  truncated  to  account  for  the  finite  mirror  limits.  Finally, 
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the  phase  of  the  field  is  reversed  to  account  for  the  change  in  direction 
caused  by  the  beam's  reflection.  The  diffractive  propagation  of  the  beam 
is  performed  using  the  Huygens-Fresnel  scaler  approximation. 

The  Huygens-Fresnel  approximation  to  the  Kirckhoff  integral  involves 
the  convolution  of  a  known  kernel  with  the  input  field  to  yield  the  out¬ 
put  field  distribution.  It  should  be  noted  that  almost  the  entire  body  of 
existent  diffraction  theory  can  be  traced  back  to  this  relationship. 

Fourier  theory  is  used  to  perform  the  convolution  since  the  Fourier  trans¬ 
form  of  the  output  field  can  be  expressed  as  the  product  of  the  Fourier 
transform  of  the  input  field  multiplied  by  the  Fourier  transform  of  the 
kernel.  Within  the  Huygens  approximation,  the  transfer  function  of  the 
Fourier  transform  of  the  kernel  is  known  in  analytic  form.  The  Bell  ap¬ 
proach  proceeds  by  numerically  Fourier  transforming  the  input  field  dis¬ 
tribution  and  multiplying  by  the  transfer  function.  The  desired  output 
field  after  propagation  is  then  determined  by  an  inverse  Fourier  transform 
of  the  resultant  product. 

The  Fourier  transform  is  obtained  using  a  numerical  technique  known 
as  the  FFT.  For  a  given  input  vector  having  N  elements,  the  FFT  will  re¬ 
turn  N  output  elements  which,  under  given  circumstances,  may  be  interpreted 
as  being  an  approximation  to  the  Fourier  transform  of  the  input  vector. 
Great  care  must  be  exercised  in  using  and  interpreting  the  results  of  the 
FFT.  The  function  to  be  transformed  is  sampled  at  uniformly  spaced  loca¬ 
tions.  These  sample  points  together  with  a  number  of  zeros  form  the  input 
vector  to  the  fast  Fourier  transform.  The  total  number  of  input  vector 
points  together  with  the  ratio  of  real  to  dummy  arguments  determines  the 
returned  maximum  extent  and  spacing  of  the  transformed  function  in  fre¬ 
quency  space.  During  these  studies,  the  Bell  code  used  an  FFT  input  vec¬ 
tor  of  21'^  of  which  only  900  points  represented  the  actual  optical  field 
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to  be  transformed.  The  remaining  points  were  loaded  with  zeros.  This 
arrangement  may  be  shown  to  theoretically  allow  for  a  maximum  frequency 
space  span  of  ±900  cycles  per  mirror  diameter  with  resolution  of  approxi¬ 
mately  0.1  cycles.  It  may  also  be  shown  that  for  the  resonator  under 
study,  an  ideal  point  sampling  scheme  would  have  utilized  an  even  higher 
ratio  of  total  FFT  points  to  optical  field  points.  At  present  it  is  suf¬ 
ficient  to  note  that  this  is  a  key  factor  in  determining  the  maximum  ef¬ 
fective  spatial  frequency  range  of  the  calculat ional  method. 

The  Air  Force  optics  code  also  begins  with  an  arbitrary  field  distri¬ 
bution  and  propagates  it  around  the  resonator.  In  the  collimated  beam 
direction  (i.e.,  from  the  large  to  the  small  mirror),  the  field  is  expressed 
on  a  plane  perpendicular  to  the  optical  axis.  At  the  small  mirror,  a 
change  is  made  to  (r,cf>)  coordinates  with  the  radius  chosen  to  match  the 
curvature  of  the  expanding  geometric  mode.  By  assuming  that  the  angular 
and  cartesian  points  are  approximately  coincident,  the  curvature  of  the 
mirror  is  absorbed  into  the  coordinate  change.  Additional  correction  need 
only  be  made  to  account  for  the  finite  absorption  and  dimensions  of  the 
mirror  together  with  the  phase  change  due  to  reflection. 

The  conversion  to  this  expanding  coordinate  scheme  with  its  associated 
elimination  of  the  multiplicative  mirror  curvature  term  is  advantageous 
since  the  spatial  frequency  content  of  the  optical  field  is  reduced,  there¬ 
by  reducing  both  real  and  frequency  space  sampling  requirements.  One  minor 
drawback  to  the  coordinate  change,  however,  occurs  because  of  the  desire  to 
keep  a  constant  number  of  sample  points  in  both  collimated  and  expanding 
space.  For  this  reason  it  is  necessary  for  the  code  to  interpolate  at  the 
small  mirror  to  provide  a  number  of  intermediate  field  samples  to  replace 
those  points  that  were  lost  in  forming  the  output  portion  of  the  beam  in 
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collimated  space.  The  interpolation  is  provided  by  fitting  the  complex 
field  function  with  a  spline  and  extracting  the  required  intermediate 
data  points. 

An  analysis  of  the  general  expanding  vs  fixed  mesh  techniques  in¬ 
dicated  that  the  observed  output  variations  could  not  be  explained  merely 
on  the  basis  of  the  differences  in  the  coordinate  system  representations 
or  the  effects  of  the  associated  interpolation,  assuming  both  codes  are 
functioning  properly.  In  order  to  remove  differences  caused  by  FFT  sam¬ 
pling  variations,  an  additional  propagator  was  written  and  added  to  the 
Bell  code.  The  test  propagator  solved  the  basic  Fresnel -Kirckhoff  in¬ 
tegral  using  standard  Fresnel  integral  subroutines  with  the  convolution 
performed  numerically.  This  basic  approach  was  formulated  in  the  conven¬ 
tional  form  for  planer  propagation  which  is  suitable  for  operation  within 
the  framework  of  the  Bell  code.  Care  was  taken  that  sufficient  sample 
points  were  available  for  an  accurate  evaluation.  The  results  agreed  with 
the  output  obtained  using  the  original  Fourier  transform  propagator  in  the 
Bell  code. 

The  numerical  aspects  of  the  particular  case  under  study  were  examined. 
It  was  hypothesized  that  variations  in  the  spatial  frequency  content  would 
vary  as  the  iterative  process  proceeded  causing  possible  numerical  problems 
at  intermediate  steps.  This  effect  was  treated  both  theoretically  in  terms 
of  determining  the  required  number  of  points  and  experimentally  by  varying 
both  the  number  of  points  and  the  FFT  guard  band.  Both  of  these  approaches 
indicated  that  the  observed  effect  was  not  associated  with  propagator  nu¬ 
merics.  A  similar  analysis  indicated  sufficient  field  sampling  points  such 
that  the  additional  mirror  curvature  term  in  the  Bell  model  was  well  re¬ 
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After  eliminating  these  potential  sources  for  the  observed  discrepan¬ 
cies,  the  AFOPT  averaged  kernel  propagator  was  considered  as  the  most  likely 
culprit.  In  order  to  check  this  portion  of  the  code,  the  AFOPT  code  was 
temporarily  modified.  The  averaged  kernel  in  the  AFOPT  code  does  not  have 
a  known  analytic  transfer  function  and  therefore  is  numerically  determined 
using  an  FFT.  This  numerical  portion  of  the  code  was  bypassed  and  replaced 
by  an  equation  representing  the  analytic  form  of  the  unaveraged  transfer 
function  in  both  the  collimated  and  expanding  cylindrical  coordinates.  In 
addition  to  added  efficiency,  this  approach  prevented  the  generation  of  high 
frequency  numerical  garbage  which  would  have  occurred  if  attempts  had  been 
made  to  numerically  compute  the  transfer  function  from  the  unaveraged  kernel. 

A  comparison  of  the  averaged  and  unaveraged  transfer  functions  showed 
some  slight  variation  in  the  significant  low  frequency  portions  of  the  spec¬ 
trum.  The  high  frequency  components  were  of  course  damped  by  the  averaging 
process  but  this  is  of  no  significance.  The  averaged  transfer  function 
showed  a  slight  premature  dip  just  above  zero  frequency  but  then  returned 
and  tracked  well  with  the  unaveraged  version.  Unfortunately,  this  pertur¬ 
bation  was  found  to  have  a  negligible  effect  on  the  calculated  intermediate 
iterative  results  in  the  AFOPT  code. 

The  final  source  of  the  variation  considered  was  the  AFOPT  spline  in¬ 
terpolation  scheme.  When  interpolation  takes  place  at  the  small  mirror, 
additional  points  are  loaded  with  interpolated  values  to  compensate  for 
those  lost  due  to  the  continuously  expanding  coordinate  system.  It  was  con¬ 
jectured  that  this  finer  resolution  together  with  its  smoothed  interpolated 
data  fit  might  tend  to  surpress  high  frequency  components.  This  was  tested 
by  overriding  the  interpolation  routine  and  filling  intermediate  grid  points 
with  neighboring  values.  By  this  technique  the  field  resolution  would  be 
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effectively  the  same  as  that  used  in  the  Bell  code.  This  again  proved  to 
be  a  null  experiment  in  that  essentially  no  change  in  the  output  was  ob¬ 
served  . 

At  the  conclusion  of  the  contract,  the  source  of  the  discrepancy  had 
not  yet  been  isolated  although  a  number  of  possible  causes  had  been  elimi¬ 
nated.  If  time  had  permitted,  the  next  step  would  have  involved  a  step¬ 
wise  investigation  of  the  actual  initialization  process  used  in  the  codes. 

Since  the  BATOPT  and  AT'OPT  codes  gave  the  same  empty  cavity  result  at 
convergence  even  though  the  intermediate  results  differ,  several  of  our 
previous  loaded  cavity  cases  were  run  to  convergence  with  the  AFOPT  + 

NORO-I  code  to  determine  the  role  of  the  gain  medium.  The  single  line 
J  =  8  case  for  a  confocal  unstable  resonatoi  with  50»  geometric  outcoupling 
was  run.  The  results  are  compared  to  our  previous  results  for  this  case 
which  were  obtained  using  the  Bell  Aerospace  Textron  optics  code.  Figs.  38  - 
42.  As  can  be  seen,  the  results  from  the  two  codes  are  in  excellent  agree¬ 
ment  with  each  other.  Since  the  saturated  gain  fills  the  resonator,  a  con¬ 
verged  solution  is  reached  in  about  IS  iterations.  For  a  second  check  case, 
the  single  line  J  =  4  case  for  the  same  resonator  was  run.  In  this  case, 
the  saturated  gain  only  fills  half  the  resonator  which  results  in  an  oscilla¬ 
tion^.  The  convergence  parameters  versus  iteration  number  are  shown  in 
Figs.  43  and  44.  The  AFOPT  +  NORO-I  code  gives  a  6-7  iteration  oscillation 
whereas  the  BATOPT  +  NORO-I  gives  a  3-4  iteration  oscillation.  A  check  of 
the  location  of  the  peak  of  the  average  intensity,  Fig.  45  shows 
that  while  the  intensity  peak  for  the  BATOPT  +  NORO-I  is  at  about  0.1,  the 
intensity  peak  for  the  AFOPT  +  NORO-I  is  about  0.15;  for  these  locations  of 
the  intensity  peaks,  the  4  and  6  iteration  oscillations  are  what  should 
occur.  These  differences  could  be  due  to  one  of  two  factors,  the  size  of  the 
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Figure  39.  Convergence  parameters  versus  iteration  number  for  single  line  P?(8)  for  AFOPT  +  NORO  I 
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gure  45.  Comparison  of  the  average  intensity  sent  to  the 

kinetics  given  by  BATOPT  +  NORO  I  and  AFOPT  +  NORO  1 
for  single  line  P_(4). 
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guard  band  used  (8196  for  BATOPT  and  2048  for  AFOPT)  or  to  the  different 
treatments  of  the  kernels  (average  kernel  in  AFOPT  and  exact  kernel  in 
BATOPT),  although  based  on  the  empty  cavity  checks  of  these  factors,  no 
changes  were  expected. 

To  check  the  effect  of  the  guard  band,  the  single  line  J  =  4  case  was 
run  with  8196  points  total  with  the  same  912  points  on  the  mirror.  The 
same  results  were  obtained  as  with  2048  points  total  with  912  on  the  mirror. 
Thus,  the  size  of  the  guard  band  does  not  seem  to  be  responsible  for  the 
slight  difference  between  the  AFOPT  and  BATOPT  results  for  this  case. 

To  check  the  effect  of  the  averaged  kernel,  the  single  line  J  =  4  case 
was  rerun  using  the  modified  AFOPT  which  employed  the  unaveraged  kernel. 

The  results,  shown  in  Fig.  46,  are  identical  to  those  obtained  with  the 
averaged  kernel.  Thus,  at  the  conclusion  of  the  grant,  the  source  of  the 
discrepency  between  the  two  optics  codes  for  the  case  of  the  saturated  gain 
filling  half  the  resonator  had  not  been  isolated.  As  in  the  empty  cavity 
case,  the  next  step  would  have  involved  a  stepwise  investigation  of  the 
actual  initialization  process  used  in  the  codes.  However,  it  should  be 
pointed  out  that  the  qualitative  results  from  the  two  optics  codes  agree  in 
all  cases  studied. 

To  summarize,  the  BATOPT  and  AFOPT  optics  codes  give  identical  results 
at  convergence  for  both  empty  and  loaded  resonator  cases.  They  differ  at 
intermediate  iterates  before  convergence  and,  in  the  case  of  a  loaded  resona¬ 
tor  in  which  the  saturated  gain  does  not  fill  the  resonator,  they  give 
slightly  different  locations  of  the  intensity  peak  and  thus,  slightly  different 
periods  of  the  oscillation  of  the  intensity.  These  results  and  the  extensive 
tests  conducted  to  determine  the  source  of  these  differences  demonstrated 
that  the  various  theoretical  approaches  to  the  physical  optics  model  give 
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Figure  46.  Comparison  of  the  average  intensity  sent  to  the  kinetics 
given  by  AFOPT  +  NORO  1  using  the  averaged  kernel  and  the 
exact  kernel  in  the  propagation  routine.  Single  line 
P^^),  50%  outcoupled,  confocal  unstable  resonator,  1024 

points  total  with  912  points  on  the  mirror. 
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surprisingly  good  agreement.  The  studies  with  the  spline  interpolation 
routine  appeared  to  indicate  that  a  simple  duplication  of  adjacent  values 
provided  essentially  the  same  result  as  an  elaborate  curve  fit  routine. 

In  addition,  the  simple  analytic  form  of  the  exact  transfer  function  ap¬ 
peared  to  agree  with  the  more  complex  numerical  approach  used  in  the 
averaged  kernel  technique.  These  comparisons  indicated  that  both  optics 
codes  could  be  modified  to  operate  more  efficiently. 

The  AFOPT  code  was  coupled  to  the  new  kinetics  model  developed  in 
this  study.  The  resulting  code,  AFOPT  +  MNORO,  has  been  checked  out  to 
determine  that  it  operates  correctly,  but  unfortunately,  the  grant  ter¬ 
minated  before  the  preceeding  Fabry-Perot  cases  could  be  run  with  an  un¬ 
stable  resonator. 


V.  CONCLUDING  REMARKS 


From  the  detailed  comparison  between  the  Blaze  II  and  MNORO  results, 
the  objective  of  developing  an  efficient  rotational  nonequilibrium  model 
of  a  cw  chemical  laser  which  is  capable  of  giving  quantitative  predictions 
of  the  power  spectral  distribution  has  been  achieved.  This  model  has  been 
coupled  to  the  AFWL  strip  resonator  code.  Unfortunately,  the  grant  ended 
before  the  coupled  code  could  be  exercised.  In  the  process  of  this  develop¬ 
ment,  several  problem  areas  have  come  to  light. 

In  addition  to  the  recognized  difficulties  with  the  HF  rate  package, 
the  prediction  of  lasing  on  the  upper  vibrational  bands  and  on  3  +  2  for 
cases  in  which  the  experimental  data  indicates  no  lasing  above  2  1  sug¬ 

gests  a  fundamental  problem  with  the  HF  rate  package,  either  in  the  pumping 
rates  to  the  higher  vibrational  levels  or  in  the  deactivation  rates  of  those 
higher  levels.  Until  this  question  is  answered,  the  determination  of  mixing 
parameters  by  requiring  total  power  and  lasing  zone  length  to  agree  with 
experimental  data  is  suspect. 

The  effect  of  F  atom  wall  recombination  is  to  slow  down  the  mixing. 
Does  F  atom  wall  recombination  really  occur  or  is  this  a  way  of  correcting 
errors  in  the  mixing  model? 

What  all  of  these  comments  indicate  is  the  need  for  a  series  of  ex¬ 
periments  to  determine  the  role  of  the  kinetics  and  the  mixing.  It  is  es¬ 
sential  that  these  experiments  be  carried  out  for  a  geometry  for  which  the 
fluid  dynamic  mixing  can  be  accurately  modeled,  i.e.,  laminar  two  stream 
mixing.  For  this  geometry,  sufficiently  detailed  fluid  dynamic  measurements 
must  be  made  to  verify  the  mixing  model;  for  the  same  geometry,  power,  power 


spectral  distribution,  and  intensity  distribution  on  each  line  for  a 
Fabry-Perot  cavity  should  be  measured.  Such  a  data  base  would  permit  a 
resolution  of  the  questions  which  have  been  raised. 

The  present  model,  in  which  the  fluid  dynamic  variables  are  input 
quantities,  generalized  to  include  three  lasing  levels  would  be  quite  useful 
in  checking  the  kinetics  model  against  srch  data. 

The  detailed  comparison  of  the  BAT  and  AFWL  strip  resonator  codes 
has  provided  confidence  in  both  optics  models  and  has  suggested  ways  the 
efficiency  of  both  codes  can  be  increased.  The  coupled  AFOPT  +  MNORO  model 
should  be  exercised  to  determine  the  combination  of  parameters  such  as  num¬ 
ber  of  points  the  cavity  is  divided  into  and  the  kinetics  integrator  error 
parameter  that  minimizes  run  time  without  loss  of  accuracy. 
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